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Abstract 

We consider stationary consensus protocols for networks of dynamic agents with switching topologies. 
The measure of the neighbors' state is affected by Unknown But Bounded disturbances. Here the main 
contribution is the formulation and solution of what we call the e-consensus problem, where the states 
are required to converge in a tube of ray e asymptotically or in finite time. 

1 Introduction. 

Consensus protocols are distributed control policies based on neighbors' state feedback that allow the coordination of 
multi-agent systems. According to the usual meaning of consensus, the system state must converge to an equilibrium 
point with all equal components in finite time or asymptotically [TJ [4l [6j [9l 1101 1111 1121 1131 115] . 

The novelty of our approach is in the presence of Unknown But Bounded (UBB) disturbances [2] in the neighbors' 
state feedback. Actually, despite the literature on consensus is now becoming extensive, only few approaches have 
considered a disturbance affecting the measurements. In our approach we have assumed an UBB noise, because it 
requires the least amount of a-priori knowledge on the disturbance. Only the knowledge of a bound on the realization 
is assumed, and no statistical properties need to be satisfied. Moreover, we recall that starting from [2], the UBB 
framework has been used in many different fields and applications, such as, mobile robotics, vision, multi-inventory, 
data-fusion and UAV's and in estimation, filtering, identification and robust control theory. 

Because of the presence of UBB disturbances convergence to equilibria with all equal components is, in general, 
not possible. The main contribution is then the introduction and solution of the e-consensus problem, where the 
states converge in a tube of ray e asymptotically or in finite time. In solving the e-consensus problem we focus on 
linear protocols and present a rule for estimating the average from a compact set of candidate points, say it lazy rule, 
such that the optimal estimate for the ith agent is the one which minimizes the distance from Xi. 

The system under consideration consists of n dynamic agents that reach consensus on a group decision value by 
implementing distributed and stationary control policies based on disturbed neighbors' state feedback. Neighborhood 
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relations are defined by the existence of communication links between nearby agents. Here, we assume that the set of 
communication links are bidirectional and define a time-varying connected communication network. 

The presentation of the results is organized as follows. We first solve the e-consensus asymptotically and in finite 
time for networks with fixed topology (we look at them as switched systems with dwell time of length infinite). To be 
more precise, for a given protocol, we find a tube of minimum radius that the agents reach asymptotically. Trivially, 
any tube of radius strictly greater than the minimum one can be reached in finite time. We do this by introducing 
polyhedra of equilibrium points and studying their stability. The above result means that, in general, the value of 
e cannot be chosen arbitrarily small. We point out its relation with the amplitude £ of the disturbances. We also 
consider additional assumptions on the disturbance realization, beyond its inclusion in D, and show that different 
type of disturbances lead to different values for the minimum radius. For certain disturbance realizations the agents 
are shown to asymptotically reach O-consensus. The last part of this paper, extends the above results to the case of 
switching topology. For a given dwell time, we find a tube that can be reached in finite time. Higher dwell times 
imply the convergence to tubes of lower radius. 

The paper is organized as follows. In Section [2] we set up the new framework of switched networks under UBB 
disturbances and formulate the e-consensus problem (Problem [TJ. In Section Owe introduce the linear protocol and 
the lazy rule. In Section [4] we study networks with fixed topology. In Section [5] we extend the obtained results to 
networks with switching topology. Finally, in Section [S] we draw some conclusions. 

2 Switched networks. 

Consider a system of n dynamic agents F = {I, . . . ,n} and let £ be a finite set of possible edgesets connecting the 
agents in F. We model the interaction topology among agents through a network (graph) G CT ( t ) = {T,E a ^), with 
time variant edgeset -E CT (t) 6 £ , where a(t) is a switching function a : IR + — > I and X is the index set associated 
with the elements of £. Also, let us call switching time a time t such that cr(t~) ^ o~{t + ) and let us call switching 
interval the time interval between two consecutive switching times. In the rest of this section, to avoid pathological 
behaviors arising when the switching times have a finite accumulation point (see, e.g., the zeno behavior in [7]) and 
in accordance with [6j[l4], we make the following assumption (see, e.g., the notion of dwell time in [6] [5]). 

Assumption 1 The switching intervals have a finite minimum length r > 0. 

Henceforth r is referred to as the dwell time. We also assume that the edgesets in £ induce undirected connected 
not complete graphs on F. For each k G X, the network Gk = (r, Ek) is undirected if G Ek then (j, i) G Ek- The 
network Gk is connected if for any agent i G F there exists a path, i.e., a sequence of edges in Ek, (i, ii) ■ ■ ■ (i r ,j), 

that connects it with any other agent j G V. Finally, the network Gk is not complete if each agent i is connected 
(with one edge) only to a subset of other vertices Nik = {j '■ G Ek} called neighborhood of i. 
Each edge in the edgeset Ek means that there is communication from j to i. As (j,i) is also in the edgeset 

Ek the communication is bidirectional, namely, if agent i can receive information from agent j then also agent j can 
receive information from agent i. Also, Gk not complete means that each agent i exchanges information only with its 
neighbors. Here and in the following, 1 stands for the vector (1,1,..., 1) T . 



2.1 Unknown But Bounded disturbances. 

Let T be the set of switching times. For all i G T, consider the family of first-order dynamical systems controlled by 
a distributed and stationary control policy 

Xi = u ia(t) (xi,y (i) ) Vt>0, t<£T 

x>(t + ) = xi{t~) -iter (l) 

where y« is the information vector from the agents in N ia ^ with generic component j defined as follows, 

{yij if j G N ia(t) , 
otherwise. 

In the above equation, y^ is a disturbed measure of Xj obtained by agent i as 

yij = Xj + dij 

and dij is an UBB disturbance, i.e., — £ < dij < £ with a-priori known £ > 0. Hereafter, we denote by d — {dij, (i,j) G 
T 2 } the disturbance vector and by D the hypercube D — {d : — £ < dij < £, V(i, j) G T 2 } of the possible disturbance 
vectors. We assume that any disturbance realization {d(t) G D, t > 0} is continuous over time. Note that both d 
and D are independent of the topology of network G (which may change over time) as they are defined on all the 
possible pairs of agents in T and not only on the links between them. The continuity hypothesis on the disturbance 
realizations can be weakened and most of our results keep holding true. However, we hold the continuity assumption 
to make the proofs of our results simpler and more readable. 

2.2 Problem formulation. 

Before stating the problem we need to introduce the notions of equilibrium point for a given disturbance realization 
d(t), and of e-consensus. 

Definition 1 A point x* is an equilibrium point for a given disturbance realization {d(t) G D, t > 0} if there exists 
i>0 such that u CT ( t ) (x*, y^) = 0, for all i G T, for all t>i. 

According to the usual meaning of consensus, the system state must converge to an equilibrium point x* G {7rl} 
in finite time or asymptotically. Hereafter, when we refer to points of type irl, we always understand that n may 
assume any value in M and we denote by {nl} the set {x : 3ir G R s.t. irl}. 

Because of the presence of UBB disturbances convergence to {ttI} is, in general, not possible. This motivates 
the following definition of e-consensus, describing the cases where the system state is driven in finite time within a 
bounded tube of radius e, 

T = {x G R n : \xi-Xj\ < 2e, G F} . (2) 

Definition 2 We say that a protocol it CT ( t )(.) makes the agents reach e-consensus in finite time if there exists a finite 
time i > such that the system state x(t) G T for all t > t. Furthermore, we say that a protocol w CT ( t )(.) makes the 
agents reach e-consensus asymptotically, if the system state x(t) — ► T for t — > oo. 

The above definition for e = (say it 0-consensus) coincides with the usual definition of (asymptotical) consensus. 
However, for a generic e > 0, the e-consensus in finite time does not necessarily implies the convergence of the state 
x to an equilibrium x* G T. In other words, x can be driven to T and keep on oscillating within it for the rest of the 
time. 



Problem 1 (e-consensus problem) Given the switched system (J^), determine a ( distributed stationary) protocol u^n) (•) 
that makes the agents reach e-consensus in finite time or asymptotically for any initial state x{fS). Furthermore, study 
the dependence of the tube radius e on the sets £ and D and on the dwell time t. 

In the rest of this paper we focus on linear protocols, and present a rule for estimating the average from a compact 
set of candidate points, say it lazy rule, such that the optimal estimate for the iih agent is the one which minimizes 
the distance from x%. 



3 Linear protocols and lazy rule. 



A typical consensus problem is the average consensus one, i.e., the system state converges to the average of the initial 
state. Its success derives from the fact that, in absence of disturbances, it can be simply solved by linear protocols. 
Let the linear protocol be given as 

Ui<r(t)(xi,y il) ) = (Vij - Xi), for all i G V (3) 

where jjij is the estimate of state Xj on the part of agent i. For a given disturbed measure yij the state Xj and 
consequently its estimate yij must belong to the interval 

va e [yij -£,va +£]■ (4) 

The crucial point is how to select y^ from the interval [yj — f, y% + £]. The next example shows that there may 
not exist equilibria if we choose simply y^ = j/y. 

Example 1 A three-agent network with a fixed edgeset E^, N 1 k = {1, 2}, A^fc = {1, 2, 3} and = {2, 3}. A simple 
criterion is to let yij ~ yij = Xi + dij 

±1 = (x 2 + d 12 ) — Xi 

x 2 = [{xi + d 21 ) - x 2 ] + [(x 3 + d 2S ) - x 2 ] 
X3 = {x 2 + d S2 ) — X-i 

Find equilibria by imposing x = and obtain 

xi = x 2 + di2 (5) 
= d 12 + d 21 + d 23 + d 32 (6) 
x 3 = x 2 + d 32 (7) 

There exist equilibria only if d\ 2 + d 2 \ + d 23 + d 32 — 0, that is, for generic values of d± 2 , d 2 \, d 23 , and d 32 we 
cannot guarantee the convergence of the system. 

Let y^ = {yij, i £ Ni} be defined according to the lazy rule 

y W =arg. min | J2 (fti" 3 *)!- ( 8 ) 

3£ N i<r(t) 

Note that as u(., .) in protocol ((3| depends on X^jgjv. (t) fei> the existence of multiple solutions y^ for ([8]) is 
not an issue. This is clearer if one observes that multiple solutions induce the same value X^ gJv . ( ^ yij for u(., .) 



in protocol ([3|. Given the lazy rule (|8]l. protocol (O turns out to have a feedback structure as, for each i G F, the 
quantity X^jgjv. V^i can be computed as 



E j6JV ^ (t) VH + \ N i*(t)\£ if 



+£ 



\N ia ( t )\xi 



- £ < Xi < 



+ £ 



(9) 



T,jeN ic . M - \ N ^(t)\£ if Xi < ^'i^'^'i - i 



Hereafter, when we refer to the linear protocol ©, we always understand that the agents choose y^ as in 



4 Fixed topology. 

In this section we consider a network with fixed topology, i.e., a network G = (T,E), with edgeset E constant over 
time. As the edgeset E remains constant, for the easy of notation, we drop the index a(t) from all the notation 
used throughout this section. Also when we refer to system JT} and to a protocol $Jfy we always mean that they are 
associated to the network G. 



4.1 Equilibrium points. 

For a network with fixed topology, we prove that the equilibrium points exist and belong to polyhedra depending on 
the type of disturbance realization. In particular, we state a first result in the case of constant disturbance d, and 
extend such a result to the case where the disturbance d takes on values in specific subsets of D. 

Lemma 1 Given the system (U) on G = (F,E), implement a distributed and stationary protocol u(.) whose compo- 
nents have the feedback form If the disturbance d is constant over time, then: 

(i) a point x is an equilibrium point for u(.) if and only if it belongs to the polyhedron 

P(d,E) = j*. —-t<—p^-- Xi < __ + £, Vl erj, (io) 

(ii) P(d,E) includes all the points in {nl}; in addition, {nl} = C\ deD P(d, E). 
(Hi) P(d,E) has 1 as only extreme ray up to multiplication by a non-zero scalar. 

Proof. (i) A point x is an equilibrium point if and only if m{xi,y^) = for all i G Y. This condition is equivalent 
to mingy e[Vij _z tyij+e]t jeN . | Ejgjv.fe ~ a*) I = that, as yij = xj + dij, in turn becomes 

52&j-xi) = o, Vier (ii) 

Xj + dij-Z < Va < Xj + dij + £, Vi G T, Vj G Ni. (12) 

The polyhedron P(d,E) is the projection of the solutions (xi, j/' 1 '), for i G T, of system (|ll|l - (|12| l in the space of the 
x variables. 

(ii) For any x G {ttI} it holds that Sj£ ^j Xj - Xi = 0. Also, - ^f^i ~ ~£ < ® < -^f^p 1 + £ because 
-£ < d i3 < £ for any i G T, j G Ni. Then, {ttI} C P(d,E). To prove that {ttI} = f] deD P(d,E) we show that 
P(£,E) = {7rl}. To see this last argument, from (JTUJ with dij = £ for all i and j we have that J ^ 3 — Xi < 



for all i G T and for any x G P(£,E). The latter means that the state Xi of each agent i must not be less than the 
average state of its neighbors in Ni and this situation occurs only if all the agents have the same state. 

(Hi) The vector 1 is an extreme ray as it is immediate to verify that if a; £ P{d,E) then x + nl G P(d,E) for 
any ir 6 R. To prove that a vector 1 is the unique extreme ray, up to multiplication by a non-zero scalar, consider a 
vector v not parallel to 1. We note that G P(d, E) and we prove that for some n G R the point + nv P{d, E). 
As — — ^T^l ' 3 ~ £ an( i ~~^JN A \~~~^~ + £ are nxe d values, we have that nv G P(d,E) for any n G R if and only if 

Vi J |^ v | 3 — for all i G T. The latter conditions define a linear system with n — 1 independent conditions 

(provided that G is connected) and the solutions are of type v = rjl for t)6E contradicting the hypothesis that v is 
not parallel to 1. 

□ 

In the proof of the previous theorem, we have observed that — ^ J rffi ' J — £ < < — ^ J ' J +£ for all i G T. When 
such inequalities hold strictly, P(d, E) is a full-dimensional polyhedron. Actually, any x of type (0, . . . , 0, S, 0, . . . , 0) 
belongs to P(d, E) if we choose 8 > sufficiently small. However, not all the polyhedra P(d, E) are full-dimensional 
as it is apparent by reminding that P(£, E) = {~rrl}. 

In the following we generalize the results of Lemma [1] to the case in which the disturbance is not constant over 
time. In other words, we are concerned with the study of the equilibrium points for generic disturbance realizations 
{d(t) G D, t > 0}. First, we can say that only the points in {irl} are equilibrium points for all the possible disturbance 
realizations {d(t) G D, t > 0}. To see this, observe that i) they are the only equilibrium points if d(t) = £ for all t 
and ii) condition (|8]) implies u(xi, y^) = 0, for all i G T for any realization {d(t) G D, t> 0}. 

We will show in the next lemma that, under certain assumptions, all equilibrium points belong to P(Q, E) — 
UdeQ E), for any subset Q C D. Before introducing the lemma, consider, without loss of generality, the box 
Q — {d G D : d~ < d < d + } C D where d~ and d + are in D and d~ < d + componentwise. Then it holds 

P(Q,E) = | ae R-:-§^- e <S^_ X4 <_§^ + f sV46 r}. (13) 

To prove (|13|l . denote by S the set on the rhs of (|13p and note that it holds either S D UdgQ P{d, E) and 3 C 
\J deQ P(d, E). Actually, for any d G Q, it holds d~ < d < d + then P{d,E) C 3, hence 3 D \J d£Q P(d,E). Also, to 
prove 3 C UdgQ E), consider a generic point x G 3. It belongs to P(d, E), where for any i G T we set 



dij = < 



4 i/ieJVt. %^-i.<0 • (14) 

d^t- otherwise 



As d G Q by construction, we have P(d, E) C UdgQ ^(^i E) which implies 3 C UdgQ E). Then, we can conclude 
that H3} holds true. 

In particular, it holds 



P(D,E) = I -2g < ^fj- ' - ^ < 2g, Vi G T L (15) 



Let us define, for a given realization d(t) and a subset Q of D, the value 

(i(Q, ti,t 2 ) = max{A : ti < t < t + A < t 2 s.t. d(t) G Q for all t <t < t + A}. (16) 

In other words, given a time interval [ti,<2], the value f±(Q,t\,t2) is the length of the longest subinterval where d(t) 
remains in Q. 



Lemma 2 Given the system £2) on G = (F,E), implement a distributed and stationary protocol tt(.) whose compo- 
nents have the feedback form Consider a disturbance realization {d(t) £ D, t > 0} and box Q — {d £ D : d~ < 
d < d + } C D. Assume that there exist two nonnegative finite numbers M and 5, such that fi(Q,t,t + M) > 8, for all 
t > 0. Then, equilibrium points x exist and belong to P(Q,E). 

Proof. We first observe that the points {ttI} are equilibrium points for a given disturbance realization d(t) and 
also that they belong to P(Q,E). We then prove by contradiction that x P(Q,E) cannot be an equilibrium point. 

If x P(Q,E), at least for one of its component, say it i, it holds that either — ^fff^\ — ~ — x i < J fjv'i ' J ~ £ or 

^ J | E ^ — - — Xi > — — ' 3 + £. The previous conditions imply that the value of u(xi,y(i)) is either strictly less than 
zero or strictly greater than zero for all d £ Q. Then, for any t > 0, there exists a time interval of length greater than 
or equal to S such that u(x,y^'') is always either strictly greater than or less than 0. Hence x is not an equilibrium 
point. 

□ 

An immediate consequence of the above lemma is the following corollary. 

Corollary 1 Given the system £7p on G — (T,E), implement a distributed and stationary protocol u(.) whose 
components have the feedback form Consider a disturbance realization {d(t) 6 D, t > 0} and a finite set 

Q = {Qi,Q2, •••} of boxes of D. Assume that there exist two nonnegative finite numbers M and S, such that 
fj,(Q r ,t,t + M) > 5, for all Q r € Q, for all t > 0. Then, equilibrium points x for «(.) exist and belong to 
n QreQ P (Qr,E). 

In addition, Corollary [T]gives us a hope that if the disturbance realization enjoys some general properties the system 
can reach an equilibrium point close to the set {ttI}. As an example, consider a disturbance realization in Corollary 
[1] characterized, at least, by Q = {Qi, Q2}, with Qi = {d 6 D : — £ < d < — d} and Q2 — {d G D : d < d < £}, with 
< d < £ , we obtain that the only equilibrium points x are in 

P (Qi,E)nP(Q 2 ,E) = {^f^-e<^^p-^<-^f^ + e, Vier}. (17) 

The above set obviously defines a neighborhood of the set {wl}, as {nl} C P(Q,E) for any possible subset Q of 
D. Interesting is that the radius of the neighborhood becomes smaller and smaller as d — > £ and that P(Qi,E) fl 
P(Q2,E) = {nl} if d = £. The same results hold, for all the situations in which we can guarantee the disturbance 
realizations characterized by Q — {Qi, Q2}, such that P(Qi, E) n P(Q2, E) is equal to a neighborhood of {irl} with 
a small radius. 

The following corollary asserts that P(D, E) is the minimal set including all the possible equilibrium points for a 
policy u(.) given an unknown but bounded disturbance in D. 

Corollary 2 Given the system on G — (T,E), implement a distributed and stationary protocol u(.) whose com- 
ponents have the feedback form (3$). If the disturbance is unknown but bounded in D, then 

(i) given any x £ P(D,E) then there exists a disturbance realization {d(t) £ D, t > 0} that has x as an equilibrium 
point; 

(ii) given any disturbance realization {d(t) G D, t > 0} then all its equilibrium points belong to P(D,E). 



Proof. (i) Any generic point x G P(D, E) is trivially an equilibrium point for the corresponding realization d(t) = d, 
t > 0, where d is defined as in (I14|l , 

(ii) Any disturbance realization {d(t) G D : t < 0} has, for any finite M > 0, n(D,t,t + M) = M then all its 
equilibrium points for u(.) belong to P(D, E) by Lemma [2] 

□ 

The following example shows that the value of the parameter e defining tube T in (|31f) cannot be chosen arbitrarily 
small. 

Example 2 Consider the network G = (r, £7) T = {1, . . . , n} and E = {(i, i + 1) : i — 1, . . . , n — 1}. Lei 

= a;i(t) + £ /or any arbitrary value of xi(t). This point is an equilibrium as long as dij(t) = for all i G T 
and t > 0. /n i/iis situation, the value e defining T in {22P * s equal to ^r^C- 



Corollary [2] suggests a way to determine a strict upper bound e for e. We have 

e = max max {xi — x,}, (18) 
ijer xep(d,e) 

whose brute force computation requires the solution of n(n— 1) linear programming problems of type ma,x xeP ( DtE ){xi — 
Xj}. Then, the computation of e becomes polynomial. 

4.2 Stability. 

In this subsection we prove the asymptotic stability of the equilibrium points. To this end, we have to introduce a 
basic property of the stationary protocol u(.) whose components have the feedback form (|3}- In the following, we 
denote by sign : R — + { — 1, 0, 1} the function that returns 1 if its argument is positive, -1 if its argument is negative, 
if its argument is null. 



Lemma 3 Given the system £7]) on G = (T,E), implement a distributed and stationary protocol u(.) whose compo 



nents have the feedback form 0). Either sign(ui(xi,y^)) = signQ^ ,- e w- ( x i ~ x i)) or sign(ui(xi,y^)) = 0, for each 



i 6 T, for each t > 0. 

Proof. For each i € T, for each t > 0, given the protocol Ui[xi,y^') = Y2j eN .(y"ij — x i), consider the solution of the 
linear problem that defines the value of y^ 

Zi = _ win Kr \y2(Vij-Xi)\. (19) 

If Zi — the lemma is proved. If Zi > two situations can occur, the value of Ejgjv- iSn ~ x i) is either strictly positive 
or strictly negative, for any jjij G [ya - + £], j G iV». We claim that if S jSJV .(j/ij — Xj) > then Ui(xi,y M ) > 
and EjeJVife - ^0 > °- If Ejejv^^J - as<) > for any y tj G [j/y - £,j/; 3 + £], j G N, then, by definition, 
Ui(xi,y^') > 0, as the chosen must belong to [j/^ — f, + £]• In addition, we have Ejgjv- ( x j + d-ij — £ — £») > 0, 
hence Ejeiv^j ~ x i) > E jG iv, (£ ~ > °> as ~Z < d ^ < £• 

A symmetric argument holds if Ej 6 jv 4 fe ~ x *) < °> for an y fei e [s/*i ~~ + £]> 3 G N i- 

□ 



Theorem 1 Given the system £7p on G = (T,E), implement a distributed and stationary protocol u(.) whose compo- 
nents have the feedback form Then, the system trajectory converges to equilibrium points in P(D,E). 



Proof. We prove the convergence to equilibrium points in P(D, E) by introducing a candidate Lyapunov function 
V(x) = | j) e E( x i ~ x ') 2 ■ Trivially, V(x) = if and only if a; £ {7i~l}; V{x) > for all x £ {wl}. We now prove 
that V(x) < for all x P(D, E). On this purpose, for V(x) we can write 



53 sign sign ( 2 
ier VjSiVi 





5Z ~ Xi ) 




j'SATj 



(20) 



From Lemma [3] if X^jgjv ( x i ~ Xi ) ~ then u%(x,y^') = 0. This in turns implies that V(x) is null if and only 
if Uk(x) — 0. The latter observation is sufficient to prove that i) the state trajectory converges to P(D,E) and that 
ii) the convergence is to an equilibrium point. Indeed, for t —* oo, we have V — > 0. Then, u — * and consequently 
x-*0. 

□ 

Theorem 2 Given the system £Ip on G — (T,E), implement a distributed and stationary protocol u(.) whose compo- 
nents have the feedback form Consider a disturbance realization {d(t) £ D, t > 0} and a box Q = {d £ D : d~ < 
d < d + } C D. Assume that there exist two nonnegative finite numbers M and 5, such that fi(Q,t,t + M) > S, for all 
t > 0. Then, the system trajectory converges to equilibrium points in P(Q,E). 

Proof. We prove the convergence to equilibrium points in P(Q, E) following the same argument used in the proof 
of Theorem [1] We only note that now we have more information on the disturbance. In particular we know that, for 
every time interval of length M, it spends at least a time 8 assuming values in Q. The explicit dependence of the 
disturbance on time makes the Lyapunov function time-varying. 

To prove the system stability we make use of the results in [8], We define the function p : R — * [0, oo) 

f 1 if d(t) e Q 
p(t) = { U (21) 

I otherwise 

It is immediate to verify that p(.) satisfies the conditions in Remark 3 in [8], in particular, there exists three finite 
values p, M, 8 > such that < p(t) < p, fc +M p(s)ds > 8 for all t > 0. We also define the function W : K" -» [0, oo) 

f if x e P(Q,E) 

W(x)={ , . J W ' (22) 

[ ™n s W £ [ IE+ d-e,x+d+i], deQ (lEjejv.fe ~ x *)\ jE.e^fe ~ x ')\\ »/ x£P(Q,E) 

Observe that W(x) = for x £ P(Q,E), whereas < W(x) < \m\ \j2j e Nt( x i _ x *)\ for a11 x & P (Q> E ) and for a11 
t > 0. Hence V^(a;) < — p(t)H /r (a;) < for all a; and all t > and, in particular, V(a;) = for all t > 0, only for 
a; £ P(Q,E). The system trajectory converges to P(Q,E). Finally, we note that the system trajectory converges to 
an equilibrium point as V(x) is null if and only if Uk{x) = 0. 

□ 

An immediate consequence of the above theorem is the following corollary. 

Corollary 3 Given the system £7p on G — (T,E), implement a distributed and stationary protocol u(.) whose 
components have the feedback form Consider a disturbance realization {d(t) £ D, t > 0} and a finite set 

Q = {Qi, Q2, ■ ■ •} of boxes of D. Assume that there exist two nonnegative finite numbers M and 8, such that 
n(Q r ,t, t + M) > 8, for all Q r £ Q, for all t > 0. Then, the system trajectory converges to f]g gS P{Qr, E). 



Finally, we can conclude that for a disturbance realization that in Corollary [T] and in Corollary [3] is characterized, 
at least, by Q = {Qi,Q 2 }, with Qi = {d G D : -£ < d < -d} and Q 2 = {d G D : -d < d < £}, with < d < 
the system trajectory converges to a neighborhood of the set {nl} with the ray of the neighborhood that becomes 
smaller and smaller as d — > £ and P(Qi, E) n P(C>2, E) = {7rl} if d = £. 

4.3 Bounds for 

In this subsection, we determine bounds for the minimum and maximum value that the components of x(t) assume 
over the time depending on the initial state a;(0) and for any disturbance realization {d(t) G D : t > 0}. In particular, 
we prove that, when we apply the lazy rule (|8}, we always obtain 

a(a;(0)) < lim x z (t) < f3(x(0)), for all i G T. (23) 

t — >oo 

As a further result, we also show that the difference between the maximum and the minimum agent states may not 
increase over the time. 

Given a network G and an initial state x(0), the main idea is to replace G by a much simpler network H composed 
by only two agents and such that the initial state of H is equal to the two maximal values of the initial state of 
network G. The result is that the maximal value assumed by the states of G is always bounded by the values assumed 
by the states of H . 

Let us denote by H = ({a, b}, {(a, fe)}) a system with only two connected agents. Let x H (t) be the state of H, 
namely, x H (t) — {x^(t), x^ (t)}. Let the components of x H (t) be subject to a constant disturbance d (t) = £. Let 
us also define ii(t) = argma,Xj e r{xj(t)} and i2{t) = argmax. : j e Y\{i 1 (t)}{xj(t)}, for all t > 0. Actually, ii(t) and i2(t) 
are the two agents with the first two maximal states. Obviously, ii(t) and i2(t) depend on time t. Analogously, define 
i„(t) = argmmj er {xj(t)} and by i„_i(t) = arg min jsr \ {in ( t )}{xj (t)}, for all t > 0. 

Lemma 4 Given the system {7]) on G = (T, E) with initial state x{0), implement a distributed and stationary protocol 
u(.) whose components have the feedback form (3$). Consider the system H — ({a, b}, {(a, b)}) with an initial state 
%a (0) — x i!(o)(0) and x^ (0) = rr i2 (o)(0). Then, for all t > 0, x^ (t) > Xiit), for all i G T and x^ it) > xt(t), for all 

»er\ {»!(«)}. 

Proof. Observe that Xi(t) is a differentiable variable for all i G F. The same property holds for Xa(t) and x^(t). In 
addition it holds that x^(t) > x^ (t) for any t > 0. 

At time t = the thesis holds by definition of values x^ (0) and x^ (0). By contradiction, assume that at some 
time instant i > the thesis is false, i.e., there exist some i,j G T such that either x^ (t) < Xi(t) or Xa(t) > Xi(t) but 
x b{^) < x j(t)' By continuity, there must also exists < t < i where one of the following conditions holds 

i) x"{t) = Xi(t), x"(t) = Xj(t), x"(t) > x k {t), for all k G Y \ {i,j}, and either x" (t + dt) < x z (t + dt) or 
x" (t + dt) < Xj(t + dt); 

ii) x"(t) = Xi(t), zf (t) > x k (t), for all k G T \ {i}, and x" (t + dt) < Xi(t + dt); 

hi) x" (t) > Xi(t), x"(t) = Xj(t), x"(t) > x k (t), for all k G T \ {i, j}, and x" (t + dt) < x 3 (t + dt); 

Consider case i). It holds x^ (t) — g^bii) ™ x a(t) < and Xi(t) = ~}2 reN .{yir(t) — Xi(t)) < and, in particular, 
yir(t) — Xi(t) < 0, for all r G Nt. As by hypothesis yf(,(t) — x^ (t) > y~i r (t) — Xi(t) for any r G T, we have x^ (t) > 
ii(t), hence the inequality x^ (t + dt) < Xi(t + dt) is false. It also holds x^ (t) = y^ a (t) — xff (t) > and Xj(t) — 



Ergjv,-^) (&>(*) ~ ^C*)) ^ Er-eJV 3VW \{i}&>(*) " ^C*)) + (%•<(*) ~ M*))- As b y hypothesis j/ 6a (i) - a; b (t) > 
—Xj(t) and ^ rgJv . t \{j}(2/j>(i) — W) < Oj we have xf (£) > Xj(f), hence the inequality xf (t + dt) < Xj(t + dt) 
is false. Hence Conditions i) cannot hold. 

We can use the first part of the above argument to prove that Conditions ii) cannot hold, and use the second part 
to prove that Conditions iii) cannot hold. 

□ 

The following corollary holds. 

Corollary 4 Given the system £7p on G = (T,E) with initial state x(0), implement a distributed and stationary 
protocol u(.) whose components have the feedback form Then: 

(i) The values assumed by the state trajectory x(t) for t —* oo satisfy the following inequalities for any disturbance 
realization d(t) 

x<„(0) < a(x(0)) < lim Xi(t) < /3(x{0)) < x n (0), for all i G T, (24) 

t— >oo 

where the bounds a(x(0)) and (3(x(0)) depend on the initial state x(0) as follows 

-wo)) = ^{^(o ^^Y^-W -e-lia -^+^w+e i (25) 

/3(x(0)) = nun {*f (0), ® +° ' ® + g + f la °*<® ~ f<® ~ * } ■ (26) 
('ii) TTie value X; x (t) is non increasing on t and the value Xi n (t) is non decreasing on t. 

Proof. To prove (i), consider the network H as denned in the proof of Lemma [4] and let us study the evolution of 
x H (t) over the time. If xf (0) — xf (0) < 2£, the system state evolves according to xf (i) = and if (i) = xf (0) — xf (i). 
If xf (0) - xf (0) > 2£, the system state evolves according to xf (t) = xf (t) + 2£ - xf (t) and xf (t) = xf (t) - xf (i), 
as long as < t < t, where t is such that xf (t) — xf (i) = 2£, i.e., t = | In x " ^ — ^ b — — -. For t > i, the system state 
evolves according to if (f) = and iff (t) = xf (t) — xf (t). Hence, for t > 0, we have 



xf (t) 



xf(0) 



xf (t) 



2 

xf(0)-xf(0)e 



g(-l+2t+ e -^) + 3: g (Q)(l- e -^) + a: f (Q)(l +e -^) 
2 



i/ xf(0)-xf(0) <2£ 

i/ xf (0) - xf (0) > 2£ and t < t 

i/ a:f (0) - xf (0) > 2£ and t > i 

if xf(0)-xf(0)<2£ 

i/ xf (0) - xf (0) > 2£ and t < i 

if xf (0) - xf (0) > 2£ and t > £ 



For f — ► oo. 



xf (<) and xf (i) converge to /3(x(0)) = min |xf (0), — 



(Q)+g£(o) 



+ £ + | In 



*q (O)-sf (0)-g 



(27) 



(28) 



Hence, 



2 i s i a *" 5 

limt-,.00 Xj(t) < /3(x(0)), for all i £ F. With an analogous argument, we can prove limt-t-oo x»(t) > a(x(0)) 

(0) f -*i„ (<>)+«♦„_! (0)+£ 



{^„(0), 



| In 1 }, 



for all i £ T. 



To prove (ii), observe that we have 



x il (0) = xf(0)>xf(t)>x n (f), 



where the first equality and the last inequality hold by definition, whereas the inequality xf (0) > xf (t) derives 
straightforwardly from the fact that xf (t) < for all t > 0. 



Corollary [4] (ii) proves that the system trajectory x(t) is bounded as t — > oo and also that the difference between 
the maximum and the minimum agent states may not increase over the time. More formally, denote by V(x(t)) = 
x ix (t) — x in it) then 

V(x(t)) > V(x(t + At)) for any t > and At > 0. (29) 

We use this last implication to introduce some additional results that will turn useful when dealing with switching 
topology systems. 

Denote by Voo = limt^co{xi 1 (t) — x< n (£)) the final value of V(x{t)). Observe that for some network G(T,E) and 
initial state x(0), there may exist some disturbance realizations {d(t) £ D : t > 0} such that even if V(x(0)) > Voo, 
the value V(x(t)) may be constant over some finite time interval before reaching its final value Voo- More specifically, 
there may exist t and At such that V(x(t)) — V(x(t + At)) > Voo- 

Example 3 Consider a networks of six agents with chain topology depicted in Fig. QJ The initial state is x(0) = 
[100, 100, 100, 0, 0, 0] T and disturbances are di2 = cfci = d-n = d$2 = t?34 = 1 and da = d^ — ^54 = dsa = des = — 1. 
Figure^ (a) shows the time plot of the evolution of the state x(t) for < t < 20 and as it can be seen, trajectories 
converge to the equilibrium x* = [63, 61, 55, 45, 39, 37] T with e = 26 (note that the initial deviation between maximum 
and minimum value of the state is 100,). Figure^ (b) displays a zoom of the trajectories for < t < 3 pointing out 
that V(x(t)) = V(x{0)) for0<t< 0.5. 

In the following we prove that we can always sample the state trajectory in such a way that the sequence of values 
for V(.) is strictly decreasing on time until the state is "almost" in P(D,E). 

To this end, given the system |T} on G = (T,E) and disturbance realizations {d(t) £ D : t > 0}, with a little 
abuse of notation, we denote by V(x) = maxigr{xi} — mmi S r{xi} the maximum difference between two components 
of vector x, for any x £ R n . Also, we denote by P(D,E) + u = {x £ R" : 3y £ P(D,E) s.t. \\x - yWac < v}, with 
v > 0, the set of points whose distance from set P(D, E) is not greater than v, according the Coo norm; and by 
V(x) = \ j)^E^ Xi ~ x i) 2 t ne Lyapunov function considered in the proof of Theorem [T] 

The following lemma holds 

Lemma 5 Given the system f7]) on G — (F,E). Let x and x in R n , ifV(x) < ^^-1) then V(x) < 7V(f), for any 
< 7 < 1. 

Proof. First we determine the bounds for the values of V(x) for x £ R n such that V(x) = V = const. Denote by 
E — {(i,j) : i < j, i,j £ F} the edgeset of the complete network induced by vertices in F. Observe that if x £ R n and 
V(x) = Vix 1 ) then 

n 2 V 2 



V{x) < max < - y ( Y : V(x) = V 



xa» 2 ^ v J/ w 8 

I (»,j)£B ) 

The last equality holds as V(x), for fixed V(x), is maximum when is maximum the number of couples of elements of 
x whose difference is equal to V(x). On the other hand, denote by E — {(i, i + 1) : i, i + 1 £ F} the edgeset of a chain 
network induced by vertices in F. Then 



mm W ( x * 
Da™ 1 2 ^— * 



V{x) > min I - \ ( Xi - xtf : V(x) = V 



2(n 



The latter equality holds as there surely exists a chain network defined by E such that Ysii j)ei?( Xi ~ x i) = V- ^he 
previous inequality holds as £V j)eE( Xi ~ x o) — ^ because E defines a connected network on G, hence there exists a 
path on E from the agent with the maximum value of the state and the agent with the minimum value of the state. 



We can now affirm that < V(x) < " \ (a) and that we have V(x) < -yV(x) if V(x) < { ^}% ■ The latter 

situation certainly occurs when V(x) < ^(ri-i) • 

□ 

Theorem 3 Given the system {!]) on G = (T,E), implement a distributed and stationary protocol u(.) whose com- 
ponents have the feedback form For each v > 0, 0<7<1 there exists a finite q(y, 7) > such that the 
values assumed by the state trajectory x(t) satisfy the following condition: either V(x(t + q(v, 7))) < jV(x(t)) or 
x(t + q(v,'y)) £ P(D, E) + v, for any t > and for any disturbance realization {d(t) £ D : t > 0}. 

Proof. We know that, when we implement a distributed and stationary protocol u(.) whose components have the 
feedback form ([3jl, the system trajectory x(t) converges to a point in P(D, E) for any x(0) £ R n . We also know that 
V(x(t)) is strictly decreasing along x(t). Then, for any x(t) £ M n , there exists a finite q(x(t),v,"/) > such that 
either V(x(t + q(x(t), v, 7))) < jV(x(t)) or x(t + q(x(t), v, 7)) £ P(D, E) + v, for any t > and for any disturbance 
realization {d(t) € D : t > 0}. As the system converges (exponentially) for any x(0) £ R n , the value q(x,v, 7) is 
finitely bounded for — > 00. Hence, the theorem is proved by defining q(y, 7) = max l6 in{ij(i, 7)}. 

□ 

The above result is strictly related to the e-consensus problem stated at the beginning of the paper (see Problem 
[1]). Actually, the above result means that the state converges in finite time to a tube of radius e = max{V(a;) : x £ 
P(D,E) + v}. 

5 Switching Topology. 

In the following, we generalize the results obtained in Section U to networks with switching topologies. Consider a 
network G CT ( t ) = (F,i? .( t )) that has a time variant edgeset -E CT (t) £ £■ We define an edgeset Et as recurrent for a 
given realization of a(t) if for all t > 0, there exists t^ > t such that cr(ifc) = k. As £ is finite there exists at least a 
recurrent edgeset for any realization of o~(t). Throughout the rest of the paper we assume that all the edegsets in £ 
are recurrent over time for any realization of o~(t), more formally 

Assumption 2 Any realization a(t) is such that, for all t > 0, there exists tk > t such that o(tk) = k for all k £ X. 

When we say that the edgeset E a ^ is time variant we understand that Assumption [2] holds and all the edegsets in 
£ are recurrent. Note that there is no loss of generality in Assumption [2] if it were false, the results of this section 
would hold for the subset £ £ £ of recurrent edgesets. 

A basic observations is that even in presence of switches V(x(t)) is not increasing on t as stated in (|29p . To see this 
note that the protocol M CT ( t )(.) induces a continuous and bounded state trajectory even for a network with switching 
topology |T]) on G a ( t ) = (r, i? CT ( t )). Then, let tk < ifc+i be two generic consecutive switching times. Corollary [4] 
applied at time t = tk instead of t = guarantees that Xi n (tk) < Xi(t) < Xi 1 (tk) for all i £ T and for all tk < t < tk+i- 
As x(t) is bounded, then u (T ( i )(.) is also bounded in the same time interval. Hence the first condition in @ implies 
that x(t) is continuous for tk < t < tk+i, whereas, the second conditions in Q imposes the continuity of the state 
trajectory in tk as x(t^) — x(t^) and in tk+i as x(t^ +l ) = x(t~^ +l ). As a consequence, for all t > 0, we also have that 
^^(O) < Xi(t) < Xi 1 (0) for all i £ F and V(x(t)) is not increasing. 



We also need to redefine the value fi(.) initially introduced as (|16p . In particular, for given realizations d(t) and 
o~(t), and a subset Q of D, we define 

a(Q, E k ,ti,t 2 ) = max{A : ti < t < t + A < t 2 s.t. d(t) G Q and o(t) = k for all t < t < t + A}. (30) 

In other words, given a time interval [£i,f2], the value u(Q,Ek,ti,t2) is the length of the longest subinterval where 
d{i) remains in Q and E a u) is equal to Ek- 

The following lemma generalizes Lemma [2] and Corollary [T] 

Lemma 6 Given the switched system £1J) on G a ( t ) = (T, E^m), implement a distributed and stationary protocol 
it CT ( t )(.) whose components have the feedback form {3]). Consider a disturbance realization {d(t) G D, t > 0} and a 
finite set Q — {Qi, Q2, • • •} of boxes of D. Assume that there exist two nonnegative finite numbers M and 5, such that 
fj,(Q r , Ek,t,t + M) > 5, for all Q r G Q, for all Ek G £ , and for all t > 0. Then the equilibrium points x exist and 
belong to fls^f Dq^q p iQr, E k ). 

Proof. Equilibrium points exists as {7rl} = (~] Ek£ £(~] d£D P(d,Ek) holds. Then, we can prove that any equilibrium 
point must belong to f) Etg£ P|q gQ P{Q r , Ek) using the argument in the proof of Lemma[2]for each couple (Q r , Ek), 
for all E k G £ and all Q r e Q. 

□ 

The results of Lemma |3] still hold in each subinterval between two consecutive switches, and then apply even in 
the switching case. 

To generalize the convergence results of Theorem [T] we need to introduce the following notations, for each Q C D: 
L(Q,Ek) = maxj, 6 p(Q iB( .) V(x) the maximum value of V(x) for points in P(Q,Ek); S(Q, Ek,v) = {x £ R" : V(x) < 
L(Q, Ek) + 2v} the set of points whose maximum difference between two components does not exceed L(Q, Ek) + 2v. 
It is worth to be noted that S(Q, Ek,v) are tubes of radius less than or equal to L(Q, Ek) + 2v and then, S(Q, Ek,v) C 
S(Q, Ey v) whenever L(Q, E k ) < L(Q, E k ). Also, observe that, by definition, it holds that P(Q, E k ) CP(<3,£ t ) + i;C 
S(Q,Ek,v). Finally, we introduce a minimum dwell time t(v,"/) — maxk£i{qk{v,"i)}- In other words, the minimum 
length of the switching intervals is equal to the maximal value over the different Ek £ £ of the times q(v, 7) introduced 
by Theorem [3] 

Theorem 4 Given the switched system {7p on G CT (t) = (T, E^m), with a minimum dwell time t(v, 7), implement a 
distributed and stationary protocol u a ^{.) whose components have the feedback form with y^ as in £3|). Assume 
that two values v > and < 7 < 1 are also given. The state is driven in finite time to the tube 

T= P| S(D,E h ,u). (31) 

Proof. We already know that the system trajectory x(t) is continuous and that V(x(t)) is not increasing. Denote by 
t s and t 3 +i with t s +i > t a + t(v, 7) two generic consecutive switching times such that a(t) — k for all t 3 < t < t s +i- 
From Theorem El we deduce that either V(x(t s+1 )) < jV(x(t s+1 )) or x(t s+1 ) G P(D,E k ) + v. The fact that V(x(t)) 
is not increasing and all the edgesets Ek £ £ are recurrent implies that there exists i > such that, if t s > t, we have 
x(t s +i) £ P{D,Ek) + v. Again, a not increasing V(x(t)) implies that x(t) g S(D, Ek,v), for all t > t s +i > t. As we 
can apply the above argument for all the edgesets Ek G £ , the theorem thesis follows. 

□ 



This last result gives an answer to the e-consensus problem stated at the beginning of this paper (see Problem [T|. 
Indeed, convergence to T, as in (|31|l . means that the agents have reached e-consensus with e = minjefc} where e*, is 
the radius of tubes S(D, Ek,v). 

Also note that the above theorem does not guarantee the convergence to an equilibrium point. Actually, switching 
systems may oscillate as shown by the following example. 

Example 4 Consider a family of chain networks G on the set of agents T — {1, 2, 3} and edgsets Ei — {(1, 2), (2, 3)} 
and E 2 = {(1, 2), (1, 3)} (see Fig. Let xi(0) = 2, z 2 (0) = 0, and x 3 (0) = 1, and {d(t) =d = const G D : t > 0}, 
with in particular di2 = dyj = d31 = 1 and d 2 i = d 2 3 = d$ 2 = — 1. If the switching time intervals are sufficiently long 
the systems trajectory oscillates in R 3 along the segment delimited by points [2,0,0] and [2,0,2]. Note that only the 
state of agent 3 changes over time. 

Finally, note that we could generalize Theorem [2] and Corollary [3] only in the assumption that, for all Q r 6 Q, 
for all Et £ £ , the values /J,(Q r , Ek,t,t + M) define sufficiently long intervals [t 3 ,t 3 +i) so that either V(x(t s +i)) is 
finitely smaller than V(x(t s +i)) of a finite value or x(t a +i) £ P(Q r , Ek) + f. In this case we have the system trajectory 
eventually assume values in C\ Ek££ Dq 6 q S(Q t , Ek, v). 

6 Conclusions. 

Despite the literature on consensus is now becoming extensive, only few approaches have considered a disturbance 
affecting the measurements. In our approach we have assumed an UBB noise in the neighbors' state feedback as it 
requires the least amount of a-priori knowledge on the disturbance. Only the knowledge of a bound on the realization 
is assumed, and no statistical properties need to be satisfied. Because of the presence of UBB disturbances convergence 
to equilibria with all equal components is, in general, not possible. Therefore, the main contribution has been the 
introduction and solution of the e-consensus problem, where the states converge in a tube of ray e asymptotically or in 
finite time. In solving the e-consensus problem we have focused on linear protocols and presented a rule for estimating 
the average from a compact set of candidate points. 

References 

[1] D. Bauso, L. Giarre, R. Pesenti, "Nonlinear Protocols for the Optimal Distributed Consensus in Networks of 
Dynamic Agents", Systems and Control Letters, vol 55, no. 11, pp. 918-928, Nov. 2006. 

[2] D. P. Bertsekas, I. Rhodes, "Recursive state estimation for a set-membership description of uncertainty", IEEE 
Trans, on Automatic Control 16 (2) (1971) 117-128. 

[3] M. S. Branicky, "Multiple Lyapunov Function and Other Analysis Tools for Switched and Hybrid Systems", 
IEEE Trans, on Automatic Control 43 (4) (1998) 475-482. 

[4] A. Fax, R. M. Murray, "Information flow and cooperative control of vehicle formations", IEEE Trans, on Auto- 
matic Control 49 (9) (2004) 1565-1476. 

[5] J. Geromel, P. Colaneri, "Stabilization of continuous-time switched systems", in Proc. of the 16th IFAC World 
Congress, Prague, Jul 2005. 



[6] A. Jadbabaie, J. Lin, A. Morse, "Coordination of Groups of mobile autonomous agents using nearest neighbor 
rules", IEEE Trans, on Automatic Control 48 (6) (2003) 988-1001. 

[7] D. Liberzon, Switching in Systems and Control, Volume in series Systems and Control: Foundations and Appli- 
cations., Birkhauser, Boston, MA, Jun 2003. 

[8] F. Mazenc, "Strict Lyapunov functions for time- varying systems", Automatica 39 (2003) 349 - 353. 

[9] L. Moreau, "Leaderless coordination via bidirectional and unidirectional time-dependent communication", in 
Proc. of the 42nd IEEE Conference on Decision and Control, Maui, Hawaii, 2003, pp. 3070-3075. 

[10] R. Olfati-Saber, R. Murray, "Consensus problems in networks of agents with switching topology and time-delays" , 
IEEE Trans, on Automatic Control 49 (9) (2004) 1520-1533. 

[11] W. Ren, R. Beard, E. M. Atkins, A survey of consensus problems in multi-agent coordination, in: Proc. of the 
American Control Conference, Portland, OR, USA, 2005, pp. 1859-1864. 

[12] W. Ren, R. Beard, "Consensus seeking in multi-agent systems under dynamically changing interaction topolo- 
gies", IEEE Trans, on Automatic Control 50 (5) (2005) 655-661. 

[13] H. G. Tanner, A. Jadbabaie, G. J. Pappas, "Stable flocking of mobile agents, part ii: Dynamic topology", in 
Proc. of the J^Zth IEEE Conference on Decision and Control, Maui, Hawaii, 2003, pp. 2016-2021. 

[14] L. Vu, D. Liberzon, "Common Lyapunov functions for families of commuting non linear systems", Systems and 
Control Letters 54 (5) (2005) 405-416. 

[15] L. Xiao, S. Boyd, "Fast linear iterations for distributed averaging", Systems and Control Letters 53 (1) (2004) 
65-78. 



Figure 1: Chain of six agents. 
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Figure 2: (a) Time plot of the state x for an initial state x(0) = [100, 100, 100, 0, 0, 0] T . Trajectories converge 
to the equilibrium x* = [63, 61, 55, 45, 39, 37] T with e = 26; (b) zoom of the first time instants (dotted 
rectangle in (a)), which highlights V(x(t)) = V(x(0)) for < t < 0.5. 
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networks of three agents. 



